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0^ . We discuss general properties of the Color Glass Condensate. We show 

that predictions for particle production in p(d)A and AA collisions derived 
from these properties are in agreement with data collected at RHIC. 
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■ In this paper we discuss the experimental signatures of the new form of 
nuclear matter - the Color Glass Condensate (CGC) in particle production 
at RHIC. Let us first see how the notion of the CGC arises in pA collisions 

^ ■ at high energy. Consider a process of inclusive particle production in pA 

qj ' collisions in a nucleus rest frame. At high energies the typical values of 

■ the Bjorken x are small. It is well-known that at small x hard processes 
develop over large 'coherence length' l c . In particular, a gluon production 
is coherent over l c ~ l/(Mx), where M is the proton's mass. For instance, 
at midrapidity at RHIC the coherence length of 2 GeV gluon is l c ~ 20 fm 
(at yfs = 200 GeV). It is much bigger than the size of the target ~ 6.5 
fm. This allows formal separation of the gluon production process into two 
parts: slow gluon emission described by the proton's light cone wave func- 
tion, and almost instantaneous interaction with the target at given impact 
parameter b described by the amplitude NQ(r,b,x), where r is the variable 
Fourier-conjugated to the gluon's transverse momentum k. In the one gluon 
exchange approximation, assuming that scattering on different nucleons is 
independent, one arrives at the formula 



X 



Ng(l, b, x) = 1 — exp 



±r 2 Q 2 s S(b)MVrv) 



(1) 



similar to the Glauber formula for the low energy hadron-nucleus scattering. 
Here Q s is a parameter with dimension of mass, S(b) is a nuclear profile 
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function and r = \r\. NG(r,b,x) can be interpreted as a forward scattering 
amplitude of a gluon dipole off a heavy nucleus. One can see from (^Q) that 
for hard gluons, such that 1/r ~ k S> Q s the scattering amplitude coincides 
with the usual perturbative expression 

N G eT \r,b,x) = r 2 7r 2 a sP T(b)xG{r,x)/(2C F ). (2) 

In the opposite limit k <C Q s the scattering amplitude becomes inde- 
pendent of its variables r, b and x as it approaches its unitarity limit. This 
is the phenomenon of saturation in nuclear and hadronic reactions [H- The 
scale Q s is called the saturation scale. Eq. (J2J) implies that Q 2 oc A 1 / 3 , 
where A is the atomic number. Thus, for a very big nucleus 4 » 1 the 
saturation scale becomes a perturbative scale Q s Aqcd which in turn 
means that a s (Q 2 ) <C 1. It can be argued that the total multiplicity of pro- 
duced gluons is dominated by gluons with the typical momentum k ~ Q s . 
Therefore, gluon production in pA collisions at high energies is calculable 
in the perturbation theory. 

Note, that the scattering amplitude JQ) was calculated in a quasi-classical 
approximation. The quasi-classical approximation corresponds to a quan- 
tum system with high occupation numbers. In terms of the QCD action 
<Sqcd this implies <Sqcd 3> 1- At small x gluons dominate over quarks. 
Therefore, we have 

\ f ol A x tr G^ (x) G^ u (x) > 1, (3) 

where the rescaled gluon field is defined as A^x) = g A®(x). On the other 
hand, a s (Q 2 ) <^ 1. Therefore, the typical gluon field of nucleus is of order 
of A^ ~ 1/g 1 . This strong gluon field at small coupling is called the 
Color Glass Condensate. This configuration is very much different from 
the perturbation theory where both the gluon (and quark) field and the 
coupling are small, and from the non-perturbative regime where both gluon 
(and quark) field and the coupling are large |JJJ- Phenomenologically, the 
CGC in a quasi-classical approximation manifests itself as a saturation of 
the scattering amplitude Ng{l,Ii,x) at small transverse momenta (!]. 

The CGC in a quasi-classical approximation can be thought of as a 
model of multiple rescatterings of a hadron in a heavy nucleus at high en- 
ergy. As such it has much in common with many other models of multiple 
rescatterings. In particular, their common prediction is the Cronin effect, i. 
e. enhancement of particle production at intermediate transverse momenta 
k in pA collisions as compared with pp scaled by the atomic number A. The 
origin of this effect is simple: a gluon traversing a heavy nucleus gains ad- 
ditional transverse momentum due to multiple rescatterings. On the other 
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hand, in a quasi-classical approximation the total number of particles is con- 
served. Therefore, if there are less particles with low transverse momentum, 
then there are more particles with high transverse momentum jSJ [3 |S1 . Of 
course, this effect is predicted to increase for heavier nuclei or more central 
collisions. 

However, a quasi-classical approximation breaks down at high energy 
since quantum evolution becomes an important process. Indeed, addi- 
tional gluon production is parametrically of order a s ln(l/x). Therefore, 
at x <C e 1 / 01 " a quasi-classical approximation is no longer valid. One might 
attempt to take the evolution into account using collinear factorization, 
which basically means incoherent production of gluons in the proton's wave 
function. As a result, proton will suffer more scatterings in a nucleus and 
the Cronin effect will increase with energy /rapidity. However, this expecta- 
tion contradicts the experimental data as we discuss later, see Fig. |21 The 
reason is that the collinear factorization scheme and the OPE break down as 
soon as multiple rescatterings are important, see e. g. This is because 
each additional scattering is a higher- twist effect. The effect of coherence 
of the parton evolution at high energies can be taken into account in the 
nonlinear evolution equations of QCD [H 1121 ITS] . These equations describe 
the high energy quantum evolution of the CGC. That is, if the scattering 
amplitude N G (r,b,x) is known at some initial value of xq, e. g. as given by 
(^Q) , the evolution equations allow calculation of the scattering amplitude at 
any x < xq. 

In the large iV c approximation the differential cross section for a gluon 
production can be written in the /c^-factorized form [3J 1141 115j 

jgg 4 _ Cj 1 Sa Sd 1 /" t2„t-72„ / v „,\ „-ifcT rr2 



d?k dy a s tt (2tt) 3 k 



J d 2 rV 2 z n G (z,Y-y)e-^ r -V 2 r N G (r,y), (4) 



where y = ln(l/x) and nc(r,b,y) is a forward gluon dipole scattering am- 
plitude off a proton. Sa and S p are cross sectional areas of the gold nucleus 
and proton correspondingly and Y is the total rapidity interval. The evolu- 
tion effects in a nucleus are enhanced by a factor of A 1 / 3 ^ 1 as compared 
to those in proton (deuteron). Therefore, na(r_,b,y) approximately satisfies 
the linear BFKL evolution equation ^B] (this is correct at not very high 
energies, when the Pomeron loops are small). The gluon dipole scattering 
amplitude can be related to the unintegrated gluon distribution function 
,x) as [Hj 

k,x) = / d 2 bd 2 re- i ^V 2 T N G (r,b,x). (5) 



a s (2vr)2 

The main property of (ft(k, x) which follows directly from the nonlinear 
evolution equations is the geometric scaling which means that 4>(k, x) = 
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(j)(k/Q s (x)), i. e. the gluon distribution becomes a function of only one vari- 
able at low x :.4 ! Here Q s (x) is the same saturation scale as in (|T|). Being 
the only dimensional parameter at low x, Q s (x) sets the scale for the gluon 
field. Eq. (J3J) implies A a ^ ~ Q s (x)/g. In course of evolution increases 
as x decreases due to increase of number of color sources. Hence, Q s {x) is 
increasing function of 1/x. It follows from the nonlinear evolution equation 
that QH : 

Cfi{x) = (^.Y A^GeV 2 . (6) 



where A ~ 0.3 ^Hl- The same gluon distribution function (f>(k, x) © enters 
expressions for the structure functions in Deep Inelastic Scattering. It allows 
to fit the initial value of x using experimental data collected at HERA. In 
®: x = 3 • 1CT 4 and A = 0.28. For RHIC and LHC it is convenient to 
write @ in the center-of-mass frame 

^=(^ V ) A e ±A ^ 1/3 GeV 2 . (7) 

The geometric scaling of the gluon distribution holds as long as the 
logarithms of energy gained in course of the BFKL evolution are bigger than 
the logarithms of transverse momentum gained in course of the DGLAP 
evolution: 

a s logQ 2 JA 2 ~ a s y > a s logk 2 /Q 2 , (8) 

which implies the geometric scaling in a wide kinematical region k < k geom = 
Qg/A ^Jj. The experimental evidences of the geometric scaling in DIS and 
heavy ion collisions are shown in Fig. ^ 

Another consequence of the quantum evolution on the unintegrated 
gluon distribution © is that its anomalous dimension 7 acquires strong 
dependence on the scaling variable k/Q s . In the perturbative regime k 3> 
&geom we get the usual leading- twist expression 4>(q,x) cx SaQI/q 2 modulo 
DGLAP corrections. However, in the saturation k < Q s , (p(q,x) cx Sa, i- e. 
7 — > 0. As we have already noted this signals the breakdown of the OPE. 
In the intermediate region Q s < k < k geom the saddle point of the BFKL 
amplitude is located at 7 1/2, which implies (f>(q,x) cx SaQs/q- Recall- 
ing, that Q s cx A 1 / 3 we find that at k < k geom the gluon distribution in a 
nucleus of atomic number A is less than A times the gluon distribution in 
a proton: 

Acf) p {k,x) Ap : U 

with /9=l/3at/c<C(5 s an d p ~ 1/6 at Q s < k < k gcom . Let me 
emphasize that the CGC takes into account two nuclear shadowing effects. 



ismd04p printed on February 2, 2008 



■5 




ZEUS BPT 97 

ZEUS BPC 95 

HI low Q 1 95 

ZEUS+H1 high Q 2 94-95 

E665 

x<0.01 



all Q 2 



p3- 
■o 



z 

73 



2- 



L 








KLN Saturation Model 

— X = 0.25 (19.6,130,21 

Energy ^ p p PHOBOS 

200 GeV UA5 • PRC 65 061 901 R (2002) 

130GeV □ Interpolation ■ PRC 65 061 901 R (2002) 

19.6 GeV A Interpolation a Preliminary 



100 



200 



300 : m \400 

N part' 



(a) (b) 

Fig. 1. (a) Geometric scaling of the total DIS cross section a(Q 2 , x) as a function 
of r = Q 2 '/Q 2 20 . (b) Geometric scaling of the total charged hadron multiplicity 
in Au-Au collisions at RHIC. Solid lines: prediction of the saturation model of 
0. Predict ion based on the collinear factorization is quite different dN / drj N paxt oc 



N, 



1/3 
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First one is a quasi-classical effect of multiple rescatterings. It necessarily 
requires higher twist effects to be included in a calculation. It predicts 
suppression at k < Q s followed by enhancement at k ~ Q s . Second one 
is a quantum evolution effect. It predicts suppression in wide kinematical 
region k < k geom both in the linear evolution region at k > Q s ('leading twist 
shadowing') and the nonlinear evolution one at k < Q s (saturation). Using 
(J1J) and (JSJ) it is easily seen that nuclear shadowing in 4>A(k, x) translates into 
suppression of particle production in deuteron-gold collisions [lUl lrJl 171151 
In Fig. [21 recent RHIC data for charged particle production at different 
rapidities and centralities |2J is shown. We see that at r\ = there is the 
Cronin enhancement of the particle production in dA as compared to pp 
in central collisions as predicted by the CGC as well as by many multiple 
rescatterings models |22j . This implies that quasi-classical approximation 
is valid. At pseudo-rapidity r/ = 3.2 corresponding to 25 times smaller x's 
than at ij = for the same transverse momentum, the evolution becomes 
essential. It manifests itself as suppression of particle production at large 
transverse momenta in p(d)A as compared to pp at higher rapidities and 
centralities. 
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Fig. 2. The nuclear modification factor: ratio of the charged hadron multiplicity in 
central (full dotes, b 3 fm) and semi-central (open dots, b 5 fm) dA collisions to 
those in peripheral ones(& rs 7 fm) rescaled by the ratio of corresponding numbers 
of participated nucleons. Lines: result of a simple CGC model of [S]. 



Not only that none of the existing conventional nuclear shadowing mod- 
els can explain both the Cronin effect and the suppression in the deuteron 
fragmentation region, but also none of those shadowing models can ex- 
plain the large value of that suppression j^Sj- in the framework of CGC 
both effects are predicted to follow from the nonlinear evolution equation 
[TU1 El ■ The value of suppression factor comes naturally as a conse- 
quence of ©• 

Since the suppression of charged particles in p(d)A and AA at forward 
rapidities at RHIC originates in a gluon shadowing © there should be 
similar suppression in open charm production. In that case the region of 
the geometric scaling, and hence of suppression, is rri[ < /% eom , where m\ = 
A: 2 + ml Another important signature of CGC is weakening of jet-jet 
correlations [53]. Indeed, since at low x a lot of particles with a typical 
momentum Q s can be produced in single nucleon-nucleon subcollision any 
two of them need not to be correlated back-to-back unless their transverse 
momenta are very large. Once all these pieces of evidence all collected 
together they will become a strong evidence for the Color Glass Condensate 
at RHIC. 

This research was supported by the U.S. Department of Energy under 
Contract No. DE-AC02-98CH10886. 
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